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Abstract 

Aladin data on fragmentation of ^^^Au projectiles are remarkably 
well reproduced by a bond percolation model. A critical behavior is 
identified on the basis of fluctuations of the largest fragment size. 

The present work is motivated by percolation studies pLj which demon- 
strated that a cumulant analysis of the largest fragment size distributions 
is a valuable tool in searching for a phase transition (critical behavior) in 
fragmenting systems. This method is applicable even to very small systems 
and therefore is well suited for applications to nuclear multifragmentation. 

The simulations were performed with a three-dimensional bond perco- 
lation model on simple cubic lattices with free boundary conditions to ac- 
count for the surface presence. Given a bond probability value p (control 
parameter) and the total number of sites Zq (system size), the probability 
distribution P{Zmax) of the largest cluster size Zmax is determined from 
a large sample of events. The statistical measures as the mean, variance, 
skewness and kurtosis contain the most significant information about the 
distribution. Of particular interest are the following cumulant ratios 
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where (Zmax) denotes the mean value, fii = {{Zmax — {Zmax)y) is the ith 
central moment, and Ki is the zth cumulant of P{Zmax)- K2 is the variance 
normalized to the squared mean, is the skewness which indicates the 
distribution asymmetry, and is the kurtosis excess measuring the degree 
of peakedness with respect to the normal distribution. In the transition 
region, these quantities obey with good accuracy finite-size scaling relations 
even for very small systems with open boundaries. This allows to identify 
universal (independent of the system size) features of Ki at the percolation 
transition. The transition in finite systems (the pseudocritical point) is 
associated with the broadest and most symmetric P{Zmax) distribution, 
which is indicated by = and the minimum value of K4, of about — 1. 
This criticality signal is approximately preserved when events are sorted by 
measurable variables correlated with the control parameter {e.g. Zbound) [I]- 




Figure 1: The cumulant ratios Ki of the P{Zmax) distribution as a function of 
Zbound- Bond percolation calculations versus the experimental data. 



The present work compares percolation predictions with the Aladin S114 
data on fragmentation of projectile spectators in ^^^Au + Au, In, Cu colli- 
sions at the incident energies of 600-1000 AMeV. Details of the experiment 



and general characteristics of the data were presented in [2]. 

Figure 1 examines the cumulants ratios Ki of the largest fragment size 
distribution P{Zmax)- The percolation results are plotted in the left dia- 
grams as a function of Zbound normalized to the system size Zq for three 
different system sizes that span over a range expected in the transition re- 
gion. Percolation events were generated for the bond probabilities uniformly 
distributed in the interval [0,1], and then sorted according to Zbound- As 
one can see, the pseudo critical point is located at Zhound/Zo = 0.84 inde- 
pendently of the system size. The experimental results are shown in the 
right diagrams for the Au + Au systems at 600 and 1000 AMeV and for the 
summed data sets (all targets and energies). Here, Ki are plotted as a func- 
tion of Zfjound- The percolation and experimental patterns of the cumulants 
are very similar. In particular, the percolation pseudocritical point is well 
reflected in the data at Zi,ound = 54. Based on this comparison, the (mean) 
system size at Zbound = 54 can be estimated as Zq = Zhound/0-84: ~ 64. 




Figure 2: Percolation predictions versus the experimental data at Zbound — 54: 
(a) the mean fragment multiplicity as a function of the fragment size (the largest 
fragment excluded), (b) the mean fragment size as a function of the fragment rank, 
(c) the multiplicity distribution of fragments with Z > 2, (d) the distribution of 
the second largest fragment charge. 

Once the system size is established, we can examine other observables 
related to fragment charge partitions. In Fig. 2 we compare percolation 
results with the data at the pseudocritical point Zf^ound = 54. Panel (a) 



shows the fragment size distribution. The model well describes the data 
over four orders of magnitude. As expected for the percolation pseudocrit- 
ical point, the distribution follows in some range the asymptotic power-law 
dependence shown by the dotted line Panel (b) shows the Zipf-type 
plot, i.e. the mean size of the largest, second largest,... r-largest fragments 
plotted against their rank r. The percolation model very well reproduces 
not only mean values but also event-to-event fluctuations. For example, the 
next panels show the multiplicity distribution of fragments with Z > 2 and 
the distribution of the second largest fragment charge. 

Similar comparisons performed at other Zbound values in the range from 
36 to 66 have also shown an almost perfect resemblance between percolation 
predictions and the data. The system sizes determined in the analysis are 
in a good agreement with the experimental estimates [3]. 

In conclusions, fluctuations of the largest fragment charge observed in 
the ^^^Au spectator fragmentations show the percolation pattern. In analogy 
to percolation, the pseudocritical point is identified at Zf,ound = 54 which 
corresponds to the He-Li isotope temperature corrected for secondary decays 
of 5.2 lb 0.4 MeV 0]. Detailed comparisons have demonstrated that the 
experimental fragment charge partitions are remarkably well reproduced by 
the bond percolation model with no free parameters and no corrections for 
secondary decays in a wide range of Z hound- In the context of the lattice 
gas model which is equivalent to bond percolation at the normal density, 
the success of percolation suggests that clusters are formed at the dense 
medium and isolated fragments are cold, in line with the "Little big bang" 
scenario of multifragmentation [5]. 
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